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producing  a  high  RCS  over  a  wide  range  of  aspect  angles . 

This  thesis  examines  the  RCS  performance  of  various 
common  retroref lector  geometries.  The  study  is  performed 
using  two  computational  electromagnetic  simulation  codes: 
a  method  of  moments  code  and  a  physical  optics  code.  The 
contour  plots  of  RCS  are  presented  for  different  geometries 
as  a  function  of  frequency.  For  retroref lectors  composed  of 
flat  plates,  the  plate  shape  is  varied  to  determine  the 
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I .       INTRODUCTION 

The  enhancement  of  the  radar  cross  section  (RCS)  of 
specific  bodies  above  their  normal  cross  section  has  several 
military  and  civilian  applications.  There  are  several 
methods  that  can  provide  this  enhancement  over  a  specific 
range  of  aspect  angles.  The  most  common  and  effective  is  the 
addition  of  radar  retroref lectors  to  the  body. 
Retroref lectors  are  simple  geometric  conducting  structures 
that  concentrate  the  reflected  wave  back  in  the  direction  of 
incidence .  They  are  capable  of  producing  a  high  RCS  over  a 
wide  range  of  aspect  angles  and  frequencies. 

Enhanced  RCS  is  desired  in  certain  commercial  and 
military  applications.  The  use  of  radar  retroref lectors  on 
small  ships  or  sail  boats  can  increase  their  RCS  so  that 
they  can  be  easily  seen  and  identified  by  other  ships' 
radars,  making  navigation  more  safe.  Another  commercial 
application  is  the  marking  of  harbor  channels  by  mounting 
radar  retroref lectors  on  buoys.  This  allows  all  ships  that 
enter  the  harbor  to  identify  and  verify  their  course  through 
the  harbor  from  safe  distances  and  depths,  especially  at 
night,  thereby  reducing  the  possibility  of  an  accident. 


In  addition,  radar  reflectors  can  be  used  for  military- 
purposes.  Perhaps  the  most  important  application  in  this 
category  is  the  decoy.  Decoys  are  objects  that  create 
pseudo- target  information  in  enemy  radar  and  thereby  confuse 
and  mislead  the  enemy.  The  payoff  is  that  a  high  cost 
weapon  will  be  expended  against  a  low  cost  decoy. 
Furthermore,  the  real  attack  can  approach  from  a  different 
direction  than  that  of  the  decoy.  Retroref lectors  are  also 
used  on  stealth  aircraft  to  increase  their  RCS  in  certain 
circumstances.  They  include  air-to-air  refueling  and  the 
masking  of  the  true  RCS  when  operating  the  aircraft  in  the 
open. 

The  purpose  of  this  thesis  is  to  examine  the  RCS 
performance  of  different  retroref lector  geometries.  The  RCS 
is  computed  using  computational  electromagnetic  (CEM)  codes, 
which  are  reliable,  low  cost  and  convenient  computer 
simulation  tools  for  scattering  analysis.  The  RCS  of 
commonly  used  retroref lector  geometries  is  calculated  at 
specific  frequencies,  using  two  CEM  codes:  (1)  the  method  of 
moments  code  Patch1  from  Sandia  Labs,  and  (2)  the  physical 


1  Computer  code  names  will  appear  in  boldface  type 


optics  code  Xpatch  from  Demaco.  These  codes  take  input  files 
with  geometry  information  and  calculate  their  RCS .  The 
results  are  loaded  into  MATLAB  for  presentation  and 
visualization.  In  many  cases  the  data  are  presented  in  two 
dimensional  contour  plots,  which  clearly  reveal  the 
retroref lector' s  ability  to  provide  relatively  high  RCS  over 
a  large  spatial  sector.  The  performance  of  different 
geometries  is  investigated  as  a  function  of  frequency.  For 
retroref lectors  composed  of  flat  plates,  the  plate  shape  is 
varied  to  determine  the  affect  of  the  plate  size  and  profile 
on  the  RCS. 

Chapter  II  provides  general  information  about 
retroref lectors :  their  types  and  general  characteristics  and 
common  applications.  Chapter  III  contains  a  description  of 
the  codes  Patch  and  Xpatch  and  information  on  the 
electromagnetic  calculations  they  perform.  Chapter  IV 
provides  the  results,  comparisons  of  data  and  comments  on 
these  results.  Finally,  Chapter  V  summarizes  this  research 
and  presents  conclusions . 


II.   RETROREFLECTORS 

In   this   chapter,   definitions   and   the   fundamental 
characteristics  of  retroref lectors  are  introduced. 
A.    RADAR  CROSS  SECTION 

1.    Definition 

The  radar  cross  section  a  of  a  target,  or  RCS  as  it  is 
commonly  referred  to,  is  defined  as  the  ratio  of  the 
target's  effective  isotropically  reflected  power  to  the 
incident  power  density.  A  general  expression  is 

Reflected  Power  /  Unit  Solid  Angle  ,    . 

a     =4tt 2— .  (1) 

Incident  Power  /  Unit  Area 

Mathematically  RCS  is  defined  by 

a  =lim  [  4:KR2^Ar   1  (2) 

where,  Es  and  E.  are  the  scattered  and  incident  electric 
field  intensities  respectively,  and  R  is  the  distance 
between  the  radar  and  the  target .  The  radar  cross  section 
has  units  of  area  but  is  frequently  presented  in  decibels 
relative  to  a  square  meter,  dBsm 

G,dBsm  =  lO\ogw(c,m2)  .  (3) 


2 .    Characteristics 

Radar  cross  section  depends  on  several  factors.  They 
include  the  frequency  and  polarization  of  the  incident  wave 
and  the  target  aspect  (that  is,  the  orientation  of  the 
target  as  viewed  from  the  radar) . 

Although  from  Equation  2  it  appears  that  RCS  depends  on 
the  square  of  the  distance  R ,     for  a  radar  in  the  far  zone 

of  the  target,  the  scattered  electric  field  intensity  Es 
varies  approximately  as  1/ R  .  Therefore,  in  the  target  far 
zone  a  becomes  range  independent. 

There  is  a  high  dependence  of  the  target  scattering  on 
the  incident  wave  frequency.  In  general,  there  are  three 
frequency  regions  where  the  RCS  of  any  body  has  distinctly 
different  features.  These  are  the  low- frequency,  resonance 
and  high-frequency  regions. 

a.  Low- frequency  Region 

In  the  low- frequency  region  the  target  size  is 
small  compared  with  the  wavelength.  That  means  if  the  target 
length  is  I  ,  then  k£<<l  where  k  is  the  wave  number,  defined 
as 
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k  ■  T  (4) 

and  X  is  the  wavelength.  In  this  region  RCS  varies  as  A.-4 
and  the  exact  shape  of  the  target  is  not  important.  This 
region  is  also  called  Rayleigh  region.  Bodies  tend  to 
scatter  isotropically  in  this  region,  although  at  quite  a 
low  level . 

b.  Resonance  Region 

In  the  resonance  region  k£~l.  For  these 
frequencies  the  induced  current  varies  in  phase  across  the 
body  and  the  RCS  is  oscillatory  with  frequency.  This  region 
is  also  called  Mie  region,  named  after  the  researcher  who 
originally  solved  the  sphere  scattering  problem. 

c.  High-frequency  Region 

In  the  high-frequency  region  the  target  size  is 
large  compared  with  the  wavelength,  that  is,  k^>>l.  Within 
this  region,  which  is  also  called  the  optical  region,  the 
scattered  field  is  strongly  dependent  on  the  angle.  Figure  1 
shows  these  three  regions  for  the  case  of  the  target  being  a 
sphere.  The  RCS  has  been  computed  using  the  Mie  series  [Ref . 
1]  • 


Figure  1.   RCS  and  frequency  regions  of  a  sphere   (After 
Ref . [1]) . 


B .    REFLECTORS 

Retroref lectors  (or  retrodirective  reflectors)  have 
been  designed  to  reflect  a  signal  towards  the  source,  over  a 
wide  range  of  incident  wave  angles .  They  are  used  in  order 
to  enhance  the  RCS  of  specific  targets  such  as  navigation 
buoys,  sail  boats,  small  drones  and  weather  balloons.  The 
important  property  is  that  they  reflect  the  incident  signal 
back  in  the  same  direction,  thus  providing  a  much  larger 
signal  than  if  the  incident  wave  were  specularly  reflected. 
An  example  is  shown  in  Figure  2 .  The  monostatic  return  from 
a  tilted  flat  plate  is  small,  but  if  a  second  plate  is 
added,  the  structure  becomes  a  retroref lector  and  the 
monostatic  RCS  is  enhanced  as  illustrated.  This  is  an 
example  of  one  of  the  most  common  types  of  retroref  lectors : 
the  corner  reflector. 

Since  the  retroref lector  reflects  the  incident  rays 
back  in  the  same  direction,  they  are  effective  only  for 
monostatic  radar  cases  and  for  special  bistatic  radar  cases 
where  the  transmitting  and  receiving  antennas  are  very 
close.   In  general,   the  monostatic  RCS  of  a  large  flat 


structure  based  on  the  physical  optics  approximation  [Ref 
1]  is  given  by 


Low 
Scatter 


High 
Scatter 


Figure  2.  The  monostatic  RCS  is  enhanced  with  the  addition 
of  a  second  plate. 


o(<t,e)«4c 


4(<M) 


(5) 


where  Ae  is  the  effective  area  of  the  reflector  presented 
to  the  source.  The  angle  9  is  measured  from  the  normal  to 
the  surface.  In  many  cases  this  formula  is  an  acceptable 
approximation  for  the  RCS  of  structures  that  are  not  flat. 

There  are  several  types  of  corner  reflectors,  but  the 
most  common  are  the  dihedral  and  the  trihedral  corner 
reflectors  which  are  discussed  in  the  following  subsections. 
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1.    Dihedral  Corner  Reflectors 

A  dihedral  corner  reflector  can  be  built  by  joining  two 
reflecting  surfaces  (such  as  rectangular  or  triangular 
plates)  at  right  angles.  An  example  of  a  dihedral  consisting 
of  two  rectangular  plates  is  shown  in  Figure  3 . 


Y 

-► 


Figure  3.  A  dihedral  corner  retroref lector  consisting  of  two 
rectangular  plates  (After  Ref .  [2] ) . 


An  incident  wave  approaching  in  a  plane  perpendicular 
to   the   intersection   line   (i.e.,  x-y      plane)   will   be 
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specularly  reflected  twice  thus  giving  a  high  reflected 
wave.  The  basic  limitation  of  a  dihedral  in  performing  as  an 
efficient  retroref lector  is  that  the  incident  wave  must  be 
restricted  near  to  the  plane  transverse  to  the  intersection 
of  the  plates  and  the  bisector  of  the  corner  angle.  That  is, 
within  a  few  degrees  of  the  x  —  y  plane  and  (J>  =  45°  in  Figure 
3. 

The  RCS  of  the  rectangular  dihedral  is  approximately- 
given  by  Equation  1  where  the  effective  area  in  the  vicinity 
of  9  =90°  is 

A=2abs'm(K/4+<b)  .  (6) 

Thus,  an  estimate  of  the  RCS  of  this  corner  reflector  is 
approximately 

\6jia2b2  sin2(7i  /  4  +(J>) 
a*    -2 .  (7) 

The  horizontal  coverage  (range  of  fy  in  a  plane  transverse 
to  the  common  edge)  of  this  type  of  reflector  is  about  90 
degrees  centered  in  the  bisecting  plane.  The  limitation  of 
the  dihedral  corner  reflectors  can  be  overcome  by 
incorporating  a  third  plate  at  right  angles  to  the  first 
two,  i.e.,  the  trihedral  corner  reflector. 
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2 .     Trihedral  Corner  Reflectors 

A  trihedral  corner  reflector  can  be  formed  by  adding  a 
third  reflecting  surface  perpendicular  to  the  first  two.  The 
spatial  sector  of  high  RCS  increases  to  a  solid  angle  of 

about  (tc/2)2  steradians  {n/2  radians  in  each  of  the  two 
principal  planes)  .  Figure  4  illustrates  three  different 
types  of  trihedral  corner  reflectors:  triangular,  circular 
and  rectangular,  along  with  approximate  formulas  to 
calculate  their  RCS .  From  these  formulas  it  is  easily  seen 
that  a  rectangular  plate  provides  a  factor  of  nine  times 
larger  RCS  than  that  of  a  triangular  plate. 

In  order  for  a  reflector  of  this  category  to  work  as  a 
retroref lector,  incident  rays  must  be  reflected  back  from 
all  three  surfaces.  To  obtain  a  larger  RCS,  the  effective 

area  At  of  the  collection  of  plates  must  be  increased.  The 
maximum  occurs  when  the  incident  wave  is  parallel  to  the 
axis  of  symmetry  of  the  reflector.  In  a  ground  based 
application  that  means  if  the  incident  wave  is  expected 
along  the  horizon  (slightly  above  the  horizontal)  ,  then  the 
trihedral  must  not  be  mounted  with  any  of  the  three  surfaces 
in  the  horizontal  plane. 


13 


Wilt 


<T  =  15.6F 


(7=3 


47tt 


Figure  4.  Triangular,  circular  and  rectangular  corner 
retroref lectors  along  with  their  formulas.  The  quantity  I 
is  the  common  edge  length  (After  Ref . [3]) . 
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The  approximate  effective  area  of  a  trihedral  can  be 
determined  by  projecting  the  surfaces  onto  a  plane 
perpendicular  to  the  viewing  direction.  For  the  case  of  a 
triangular  plate  trihedral,  Figure  5  shows  how  the  effective 
area  can  be  determined.  This  area  is  represented  by  the 
shaded  part.  After  Ae  is  determined,  the  RCS  can  be 
estimated  using  Equation  1. 


\A-''f-\^ 

"* 

"*»>*. 

<**■>•*        SJ( 

/  \'-';»3 

j^ 

^■k^    S '  *ar  \ 

<*  ?'w 

••'^•"•^P^Ov 

Figure   5.   Determination   of   the   effective   area   of   a 
triangular  plate  trihedral  (After  Ref .  [3]) . 

For  a  square  plate  trihedral,  the  RCS  can  be  calculated 
as  a  function  of  the  spherical  angles  0  and  (|)  .  For  values 
of  0  and  <|)  near  the  symmetry  axis  the  result  is 
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a(9,<j))*  ^cos2e(4-cot(t>)2r  (8) 

A, 

where  I     is  the  side  length.  In  a  similar  manner,  for  a 
triangular  corner  reflector,  the  RCS  can  be  calculated  by 

o(9, <1>)  «  ~y ^  cosG  +  sin0(sin<|)  +  cos<|))  -  2|cos0  +  sin0(sin(|)  +  cos<J))J 

(9) 

The  above  formulas  are  limited  in  application  and  even 
then  do  not  provide  accurate  results  in  all  cases.  For  more 
accurate  RCS  prediction  computer  codes  must  be  used.  Several 
codes  that  can  be  used  to  solve  this  type  of  problem  are 
discussed  in  the  next  chapter. 
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III.  COMPUTER  CODES 

A.    Patch 

Patch  is  a  method  of  moments  triangular  subdomain  code 
that  was  developed  by  Sandia  Labs.  It  is  a  Fortran  computer 
code  that  is  based  on  the  solution  of  the  electric  field 
integral  equation  (EFIE) .  It  is  applicable  to  both 
electromagnetic  radiation  and  scattering  problems .  It  models 
a  continuous  surface  by  discrete  triangular  patches  in  order 
to  handle  arbitrary  body  shapes  easily  and  accurately. 

1.    Method  of  Moments  (MM) 

The  MM  is  a  numerical  technique  that  is  used  to  reduce 
the  EFIE  to  a  set  of  linear  equations,  and  thus  converts  the 
integral  equation  to  a  matrix  problem.  The  size  of  the 
matrix,  that  is,  the  number  of  the  unknowns,  depends  on  the 
number  of  the  triangular  patches  used  to  represent  the  body. 

The  unknown  quanity  in  the  EFIE  is  the  surface  current 

distribution  J s  .     To  apply  the  MM,  J s     is  represented  by  a 
series   with  unknown  coefficients 


l=J,iJ„  do) 


n=\ 
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where  { J„ }  are  the  basis  functions.  The  selection  of  the 
basis  functions  is  important.  Pulses,  5  functions,  triangles 
and  sinusoidals  are  frequently  used  as  basis  functions.  In 

order  to  determine  the  unknown  coefficients  /„ ,  a  matrix  of 
N    linear  equations  is  formed.  By  exploiting  the  properties 

of  the  basis  functions,  the  {/„}  are  determined  and  thus  Js 
is  also  determined. 

After  the  current  coefficients  have  been  determined, 
the  radiation  patterns  and  the  scattered  fields  can  be 
computed.  If  the  MM  is  properly  applied,  the  actual  current 
series  representation  will  converge  to  the  true  current  as 
the  number  of  the  basis  functions  is  increased.  A  useful 
rule  of  thumb  is  that  the  edge  length  should  not  exceed 
0.1X  .  In  general,  the  number  of  edges  and  unknowns  is  almost 
the  same.  Small  differences  occur  due  to  the  fact  that  there 
are  more  than  two  current  coefficients  associated  with  some 
edges.  On  the  other  hand,  if  the  edge  of  a  triangle  is  also 
the  edge  of  the  body,  then  there  is  no  current  flowing 
across  it  and  therefore  no  expansion  coefficient  is 
required.   Patch   uses   overlapping   rooftop   functions   as 
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expansion  functions  J„ .  The  current  at  a  point  on  a  surface 
triangle  is  the  vector  sum  of  the  currents  crossing  the 
three  edges  weighted  by  their  distances  from  the  edge. 

Several  computer  codes  are  used  to  generate  facet 
models,  perform  the  computational  electromagnetics,  and 
present  the  output  data.  The  code  relationships  are 
summarized  in  Figure  6 .  Several  output  files  are  generated 
by  Patch.  They  include  the  RCS  as  a  function  of  9  and  4> 
which  can  be  plotted  in  MATLAB . 

2 .    Program  Build 

Build  is  a  program  that  generates  faceted  target  bodies 
by  joining  a  set  of  simple  three-dimensional  shapes.  The 
data  is  written  in  a  format  that  is  recognized  by  Patch. 
Build  is  capable  of  generating  five  basic  geometries: 
quadrilateral,  cylinder,  cone,  disk  and  sphere: 

(1)  Quadrilateral:  This  is  formed  by  specifying  the 
coordinates  of  the  four  corners  and  the  number  of  edges 
along  two  adjacent  sides. 

(2)  Cylinder:  A  cylinder  is  formed  by  specifying  four 
points.  The  first  two  points,  a  center  point  and  a  point  on 
the   circumference  define   the  top   face  and  similarly  two 
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program 
name 

input  file 

output 

name 

file  name 

*  denotes  an  arbitrary  name  subject  to 
each  code's  character  restriction 
Extensions  must  be  included  as  indicated. 


knit.x 


atp.x 


in.  facet 


*.facet 


ACAD 


\  facet 


out.  patch 


i 


tges 
k.igs 


AutoCad 


_ 


bldmatx     * 


mbuild.x 


lC  outpatch    y*~      mPatch 


i 


inpatch 


m  pattern 
files 


MATLAB 


1 


plot  patterns 


*.m  geometry 
files 


*r 


plot  geometry 
pltpatclT.i 


m 


in.  patch 


out.  facet 


pta.x 


in.  patch 


CODES 

AutoCad,  ACAD:  computer  aided 
drawing  (CAD) 

atp.x:  ACAD  to  inpatch  translator 

pta.x:  inpatch  to  ACAD  translator 

bldmatx:  build  matlab  geometry 
files 

mbuild.x:  build  simple  geometries; 
specify  excitation  parameters  for 
mpatch.x  calculations. 

mpatch.x:  MM  current  and 
scattering  calculations 

knit.x:  remove  duplicate  edges; 
check  for  geometry  file  errors 


Figure  6.  Relationships  between  computer  codes 
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other  points  define  the  bottom  face.  The  cylinder  axis  is 
defined  by  the  two  center  points. 

(3)  Cone:  A  cone  with  finite  length  and  open  end  is 
formed  by  specifying  three  points.  These  are  the  top,  the 
center  point  of  the  base  and  a  point  on  the  circumference. 
The  number  of  edges  from  the  apex  to  the  base  and  around  the 
circumference  are  also  specified. 

(4)  Disk:  A  disk  is  formed  by  specifying  a  center,  a 
circumference  and  an  axis  point.  The  number  of 
circumferential  and  radial  edges  are  also  required  in  order 
for  the  triangulation  to  be  specified.  Nonuniform  segment 
lengths  along  the  radial  dimension  can  be  specified. 

(5)  Sphere:  A  sphere  is  specified  by  a  center  point,  a 
point  on  the  equator  of  zero  degrees  longitude  and  a  point 
at  the  north  pole.  A  portion  of  a  sphere  can  be  specified  by 
providing  the  beginning  and  ending  angles  at  specific 
longitudes  and  latitudes. 

3 .    Code  Capabilities 

Apart  from  building  the  shape,  Patch  provides  some 
options  to  the  user.  Among  them  is  the  symmetry  option.  This 
option  enables  the  user  to  reduce  the  number  of  computations 


21 


and  thus  the  time  required  by  using  image  planes  placed  on 
the  x  =  0,  y  =  0  or  z  =  0  planes.  These  image  planes  can  be 
either  perfect  electric  conductors  (PEC)  or  perfect  magnetic 
conductors  (PMC) .  By  inserting  a  symmetry  plane,  for  each 
triangular  patch  an  image  one  is  created.  Other  quanities 
such  as  voltage  sources  or  distributed  impedances  follow  the 
image  properties  and  are  taken  into  account  automatically. 
Capabilities  are  summarized  in  Table  1. 


Table  1.  Summary  of  Patch  Capabilities 


Arbitrary  Shape 

Open/Closed  objects 

Modelled  by  triangular  "  patches" 

Variable  patch  density 

Front  end  for  graphical  composition 
Arbitrary  edge  multiplicity 
Non-orientable  surfaces 
Symmetry  planes  may  be  included 
Multiple  bodies 

Surface 

Basis  functions  yield  surface  currents 
Type:  Wilton  Rao 
Free  of  line  and  point  charges 
Equivalent  Thevenin  circuits 
Lumped  and  surface  impedance  loading 


Excitation 

Voltage  sorces 
Plane  waves 
Both 

Calculated  Quanities 

Surface  currents 
Far  field  patterns 
Radar  cross  section 
Field  calculations  at 
general  observation 
points 

Frequency  Domain 

Pattern  loops 
Frequency  loops 


Apart  from  the  above  options,  Patch  has  the  options  of 
entering  the  number  of   impedance   loads  and  faces  with 
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surface  impedance  loads,  printing  the  currents  and  computing 
the  field  at  an  observation  point  or  points.  All  of  these 
options  make  Patch  code  a  powerful  and  useful  tool  in  the 
RCS  calculation  of  a  variety  of  structures. 
B .    Xpatch 

Xpatch2  is  a  computational  electromagnetics  code  from 
Demaco.  It  consists  of  high-frequency  radar  signature 
prediction  codes  based  on  the  shooting  and  bouncing  ray 
(SBR)  technique  [Ref .  4]  .  In  this  technique,  a  number  of 
rays  are  transmitted  towards  a  target  and  then  traced 
according  to  geometrical  optics  (GO)  as  they  bounce  around 
within  the  target.  When  exiting  the  target,  the  scattered 
far  field  is  calculated  using  a  physical  optics  (PO) 
integration.  Therefore,  using  the  SBR  technique,  all  bounce 
contributions  are  accounted  for  in  the  PO  current 
calculation.  Several  factors  such  as  the  ray  divergence  and 
polarization  are  included  in  the  tracing.  There  is  also  an 
option  to  include  the  first-order  edge  diffraction  which  is 
calculated  by  the  physical  theory  of  diffraction  (PTD) . 


2  All  computed  data  and  information  relating  to  Xpatch  was 
provided  by  Professor  David  C.  Jenn 
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The  code  is  composed  of  three  parts:  (1)  the 
electromagnetic  part  which  has  four  individual  codes 
(xpatchl,  xpatch2 ,  xpatch3  and  xpatch4) ,  (2)  the  CAD  and 
visualization  part,  and  (3)  the  graphical  user  interface 
(GUI) .  Any  subcode  can  be  used  individually  or  as  a 
component  in  the  GUI.  For  example,  in  xpatchl  the  total 
scattered  field  is  the  coherent  sum  of  the  fields  from  the 
first  and  higher  bounce  contributions  and  the  edge 
diffraction.  The  calculation  of  the  total  field  can  be  done 
using  either  a  PO  integration  or  a  z -buffer  technique.  The 
best  choice  of  the  method  depends  on  the  target  shape .  That 
is,  if  the  target  has  a  simple  geometry  (plate,  cube, 
sphere,  etc.)  the  first  method  is  the  best  choice.  For  more 
complex  targets  (plane,  missile)  the  z-buffer  is  faster.  All 
three  methods  give  nearly  the  same  results  for  complex 
targets . 

The  codes  discussed  in  this  chapter  were  used  to 
analyze  the  RCS  of  the  different  retroref  lectors  and  the 
results  are  presented  in  the  next  chapter. 
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IV.   RESULTS  AND  DATA  ANALYSIS 

In  this  chapter  the  RCS  of  retroref lectors  comprised  of 
three  different  plate  shapes  is  computed  and  the  results  are 
discussed.  The  computation  is  based  on  both  the  Patch  and 
Xpatch  codes  and  a  comparison  of  these  two  codes  is  given. 
Furthermore,  the  code  convergence  versus  the  number  of 
segments  is  checked  and,  finally,  a  comparison  is  made  of 
the  RCS  values  from  the  codes  and  formulas  given  in  Figure 
4. 
A.    PLATE  SHAPES 

The  basic  retroref lector  shapes  that  are  considered  in 
this   thesis  are :  ( 1 )   three  rectangular  plates,   (2)   three 
triangular  plates   and   (3)   three   circular  wedge   shaped 
plates.  The  plates  are  located  in  the  x—  y,      x  —  z     and  y  —  z 

planes  as  shown  in  Figures  7,  8,  and  9.  All  of  the  principal 
edge   lengths   are  arbitrarily  set  at   10m.   All   data   is 

normalized  by  A2.  Thus,  the  RCS  is  given  by  o7A2  in  dB .  The 

A  A  A 

incident  wave  polarization  is  6     and  both  the  6      and  0 
receiver  polarizations  were  computed.  The  shapes  were  tested 


25 


10  o 


Figure  7.  Three  square  plate  geometry  with  10m  side  length 
and  2  0  edges /side. 
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Figure  8.   Three  triangular  plate  geometry  with  10m  side 
length  and  2  0  edges /side. 
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10  o 


Figure  9.  Three  circular  wedge  shaped  plate  geometry  with 
10m  side  length  and  17  edges/side. 
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at  frequencies  of  30,  60  and  120  MHz.  The  computational 
models,  that  is,  the  number  of  edges  per  side  were  generated 
using  the  CAD  program  ACAD.  The  RCS  is  computed  over  the 
hemisphere  where  z  >  0 .  This  corresponds  to  spherical  angles 
0<9<90°  and  0<(|><360o  and  direction  cosines  -l<w,v<l  where 

w  =  sin9  cos<J)  ,  (11) 

v  =  sinG  sin  (J)  .  (12) 

Contour  plots  were  generated  in  direction  cosine  space.  The 
contours  enclose  regions  where  the  RCS  is  greater  than  the 
specified  value. 

At  a  frequency  of  60  MHz,  the  RCS  is  higher  in  the  case 
of  the  rectangular  and  curved  plate  structures  (+20  dB) 
relative  to  the  case  of  the  triangles  (  +  10  dB)  .  This  was 
expected  since  the  effective  area  of  the  rectangular  plates 
is  larger  than  that  of  the  triangular  plates .  The  highest 
RCS  value  in  all  three  cases  occurs  in  the  first  quadrant  of 
the  u  —  v  plane  (0<w,v<l).  The  contour  plots  for  these  three 
cases  are  shown  in  Figures  10,  11,  and  12.  By  doubling  the 
frequency  to  120  MHz,  the  RCS  is  also  increased  for  all 
three  cases  and  the  highest  values  become  +30  and  +20  dB 
respectively.  The  area  where  the  highest  RCS  occurs   remains 
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>    0 


-0.4- 


Figure   10.   RCS   of   three  rectangular  plate   geometry  at 
frequency  of  60  MHz  and  20  edges/side  using  Patch. 
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>    0 


Figure   11.   RCS   of   three   triangular   plate   geometry   at 
frequency  of  60  MHz  and  20  edges/side  using  Patch. 
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>    0 


Figure  12.  RCS  of  three  circular  wedge  shaped  plate  geometry 
at  frequency  of  60  MHz  and  17  edges/side  using  Patch. 
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the  same  as  that  at  the  frequency  of  60  MHz.  This  was 
also  expected  since,  from  Equation  3,  the  RCS  is  inversely- 
proportional  to  the  square  of  the  wavelength. 

In  addition  to  the  three  above  shapes,  a  fourth  one  was 
also  investigated.  This  was  the  octahedral  shown  in  Figure 
13  whose  RCS  was  calculated  at  the  frequency  of  30  and  60 
MHz.  The  octahedral  can  be  considered  as  eight  trihedrals 
with  a  common  vertex.  At  the  frequency  of  60  MHz  the  highest 
RCS  value  was  +20  dB  and  this  occured  in  the  directions 
shown  in  Figure  14 . 
B.    COMPARISON  OF  RESULTS 

RCS  data  were  generated  using  both  Patch  and  Xpatch. 
This  was  done  for  verification  and  comparison  of  the 
computational  approaches.  Looking  at  the  contour  plots  from 
both  codes,  it  is  seen  that  they  generally  agree  at  high 
frequencies.  The  results  for  the  case  of  the  triangular 
plates  at  the  frequency  of  120  MHz  are  shown  in  Figures  15 
and  16.  At  lower  frequencies  the  low  level  contours  are 
different  although  the  highest  RCS  values  remain  the  same. 
For  a  detailed  look,  several  cuts  were  taken  at  angles  of 
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10  -10 


Figure  13.  Octahedral  geometry  with  10m  side  length  and  15 
edges /side . 
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Figure  14.  RCS  of  octahedral  at  frequency  of  60  MHz  and  15 
edges /side  using  Patch. 
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Figure   15.   RCS   of   three   triangular   plate   geometry   at 
frequency  of  120  MHz  and  20  edges/side  using  Patch. 
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>    0 


Figure   16.   RCS   of   three   triangular   plate   geometry   at 
frequency  of  120  MHz  and  20  edges/side  using  Xpatch. 
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0=30,  60  and  90  degrees.  A  comparison  between  the  two  codes 
for  the  triangular  plate  case  at  the  frequency  of  60  MHz  and 
for  0=60°  is  shown  in  Figures  17  and  18.  The  results  are  in 
agreement  given  that  the  two  codes  are  based  on  different 
calculation  methods   (MM   and   PO)  .   Both   the   0   and   (j> 
polarizations   from  Xpatch  are  very  close  while   the  <\> 
polarization  plot  from  Patch  is  much  different  than  that  of 
0  polarization.  This  is  due  to  the  fact  that  the  Xpatch 
calculation   did   not   include   edge   diffraction.   As   the 
frequency  decreases,  the  maximum  values  of  a  are  close,  but 
the  shapes  of  the  curves  change  dramatically.   This   is 
expected  at  low  frequencies  because  PO  and  SBR  are  high- 
frequency  techniques . 
C.    CODE  CONVERGENCE 

The  effect  of  the  number  of  facets  in  the  Xpatch  model 
was  checked.  Three  separate  cases  of  plate  segmentation  were 
examined  for  frequencies  of  30,  60  and  120  MHz:  (1)  for  25 
edges/side,  (2)  for  20  edges/side  and  (3)  for  one  edge/side 
(that  is,  three  triangles  without  segmentation) .  The  results 
for  (1)  and  (3)  are  shown   in  Figures  19   and  20.  From  the 
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Figure  17.  RCS  of  three  triangular  plate  geometry  at 
frequency  of  60  MHz  and  20  edges/side  for  0=60°  using 
Patch. 
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Figure  18.  RCS  of  three  triangular  plate  geometry  at 
frequency  of  60  MHz  and  20  edges/side  for  9  =  60°  using 
Xpatch . 
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Figure   19.   RCS   of   three   triangular   plate   geometry   at 
frequency  of  120  MHz  and  25  edges/side  using  Xpatch. 
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Figure   20.   RCS   of   three   triangular   plate   geometry   at 
frequency  of  12  0  MHz  and  1  edge /side  using  Xpatch. 
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comparison  of  the  plots  it  is  seen  that  there  are  no 
differences,  which  implies  that  the  segmentation  is  not  an 
important  factor  in  the  calculation  of  RCS  in  Xpatch  code. 
This  is  expected  because  the  "  soft"  edges  inside  of  the 
plate  do  not  contribute  to  scattering.  However,  the 
calculation  can  take  significantly  longer  (about  10  times) 
when  small  triangles  are  used  because  reflections  must  be 
computed  between  a  larger  number  of  patches.  Finally,  the 
edge  diffraction  does  not  affect  the  calculation  of  RCS  as 
shown  in  Figure  21. 

The  effect  of  the  number  of  facets  in  the  Patch  model 
was  also  checked.  Two  separate  cases  of  plate  segmentation 
were  examined  for  a  frequency  of  120  MHz:  (1)  for  25 
edges/side,  and  (2)  for  20  edges/side.  The  result  for  (1)  is 
shown  in  Figure  22 .  From  the  comparison  of  Figures  15  and  22 
it  is  seen  that  there  are  only  small  differences,  which 
implies  that  the  segmentation  is  not  an  important  factor  at 
this  frequency  in  the  calculation  of  RCS  in  Patch  code. 


43 


>    0 


-0.2 


-0.4  - 


-1    -0.8   -0.6   -0.4   -0.2 


0 
U 


0.2    0.4    0.6    0.8     1 


Figure  21.  RCS  of  three  triangular  plate  geometry  at 
frequency  of  120  MHz  and  1  edge/side  with  edge  diffraction 
using  Xpatch. 
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Figure  22.  Three  triangular  plate  geometry  at  frequency  of 
12  0  MHz  and  2  5  edges /side  using  Patch. 
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D.    COMPARISON  OF  RCS  VALUES 

The  last  step  is  to  compare  the  highest  RCS  values  that 
come  from  the  codes  with  those  from  the  formulas  of  Figure 
4.  The  results  for  the  rectangular  and  triangular  cases  are 
included  in  Table  2  and  show  that  in  general  the  values  are 
in  agreement,  especially  for  higher  frequencies.  All  of  the 


RCS  values  are  normalized  by  X2 


Table  2.  Theoretical  and  computed  maximum  RCS  in  dB . 


Plate 
Shape 

60  MHz 

120 

MHz 

240 

MHz 

codes 

formulas 

codes 

formulas 

codes 

formulas 

Rectangular 

25.54 

27.81 

39.42 

39.84 

58.17 

51.88 

Triangular 

18.55 

18.26 

31.53 

30.30 

43.87 

42.34 
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V.    CONCLUSIONS 

In  this  thesis  the  RCS  performance  of  three 
retroref lector  geometries  was  examined.  The  RCS  of  the  above 
geometries  was  computed  for  a  range  of  frequencies  using  the 
CEM  codes  Patch  and  xpatch.  The  results  showed  that  in 
general  the  two  codes  agree  in  both  the  level  of  RCS  as  well 
as  the  spatial  sectors  where  these  values  occur.  The 
agreement  increases  with  the  frequency  since  the  ray- 
approach  of  Xpatch  is  based  on  high  frequency  techniques. 
One  limitation  for  convergence  of  Patch  is  that  the  edge 
length  should  not  exceed  0.1^.  .  However,  it  was  found  for  the 
geometries  under  consideration  here  that  a  maximum  edge 
length  of  0.2>»  could  be  tolerated.  The  segmentation  is  a 
factor  that  does  not  affect  the  RCS  in  Xpatch  code. 

As  expected,  the  RCS  increases  with  the  effective  area 
of  the  retroref lector .  The  highest  RCS  values  from  these  two 
codes  are  in  agreement  with  those  from  the  theoretical 
expressions  and  this  verifies  the  utility  of  the  formulas. 
Although  there  are  some  differences  between  the  two  codes 
(given  that  they  are  based  on  different  calculation  and 
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methods) ,  the  use  of  either  Patch  or  Xpatch  is  an  accurate 
convenient  way  for  predicting  retroref lector  RCS. 

The  triangular  plate  structure  has  lower  RCS  compared 
to  the  circular  and  rectangular  plate  structures .  In  many- 
applications  the  triangular  arrangement  is  preferred  because 
of  its  smaller  volume  compared  to  the  other  two.  Smaller 
volume  has  reduced  weight,  wind  resistance,  and  is  less 
unsightly  than  a  large  package.  Thus  the  ratio 

High  RCS  Solid  Angle 
Volume 

could   serve   as   a   figure   of   merit   for   evaluating 
retroref lectors . 

The  question  arises:  how  small  can  a  retroref lector 
become  yet  still  provide  a  sufficiently  high  RCS?  From  the 
calculated  data  it  appears  that  2X  edges  are  about  the 
minimum  for  a  good  retroref  lector .  Below  that  the  RCS 
pattern  tends  to  be  relatively  isotropic.  In  the  range  of 
2X  to  4X  edge  lengths  either  of  the  simple  three  plate 
reflectors  are  efficient.  When  limited  to  this  size  range, 
other  more  complicated  shapes  do  not  significantly  improve 
the  RCS  solid  angle  to  volume  ratio. 
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